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Abstract. Some of (1 — x)NaNb03 — (x)AB03 perovskite solid solutions exhibit 
a dramatic diffusion of the dielectric permittivity e' maximum and relaxor-type be- 
havior when the second component concentration exceeds a threshold value xq. The 
concentration phase transition to the relaxor-like phase is abrupt (of the first order 
kind) that is seen from the step in the dependence of the e'(T) maximum tempera- 
ture, T m , on x. The precursor of this transition is a giant (up to 100 K) temperature 
hysteresis of s'(T). Some relaxor-like properties appear even at x < Xq in the course of 
cooling while disappear in the course of heating. The experimental data obtained are 
qualitatively described within a Landau-type phenomenological approach, assuming 
the relaxor-type behavior to be local stress-induced. 



I INTRODUCTION 

The relaxor behavior in perovskites was studied predominantly in Pb-containing 
ternary compounds (PMN, PST, PSN) and solid solutions (PLZT, PMN-PT) [1- 
3] . Last years some lead- free BaTiC^-derived solid solution compositions attracted 
much attention as environmentally benign relaxor materials [4 ]. However up to 
now the obtained values of the dielectric permittivity e' maximum temperature, 
T m , in the BaTi03-based relaxors were too low for potential applications [4]. The 
aim of the present paper is studying the dielectric properties of some NaNbC>3 - 
based solid solutions forming a new family of lead-free materials with diffuse phase 
transition (DPT). 

On heating NaNbC>3 exhibits a series of six phase transitions from the low 
temperature ferroelectric (FE) phase N (R3c) to the high temperature paraelectric 
(PE) cubic phase U (Pm3m) through different antiferroelectric (AFE) and PE 
phases [5]. An e(T) maximum originating from the first order transition between 
two AFE phases: P(Pbma) and R(Pmnm) is observed at 350-370°C. 

Similar to solid solutions of other perovskite antiferroelectrics, the NaNbC>3 - 
ABO3 binary solid solution systems can be divided into two groups [6]. In the 
solid solutions of group I (e.g. (Na,Li)Nb0 3 and (Na,K)Nb0 3 systems) the high 



temperature FE phase appears at small (a few mol.%) content x of the second 
component AB0 3 , the T m (x) dependence is rather smooth and the e'(T) maxima 
are sharp. In the solid solutions of group II the AFE phase remains stable up to a 
comparatively high x values. In contrast to the solid solutions of other perovskite 
antiferroelectrics, the T m {x) dependence of the NaNbOs-based group II solid solu- 
tions remains smooth only up to a threshold x = x value. At x > x the phase, 
usually referred to as FE (though sometimes it is supposed to be ferrielectric), 
becomes stable, which is accompanied by an abrupt drop in the T m values and 
dramatic diffusion of the e'(T) maximum. While the T m values of the composi- 
tions with x < Xq do not depend on frequency, the compositions with x > x were 
reported to exhibit a frequency dispersion of e' and an increase of T m with fre- 
quency [6]. Thus, it seems that NaNbC>3 -based solid solutions belonging to group 
II exhibit a relaxor-like behavior at x > xq- Besides, the lack of systematic data 
on the properties of such materials prevents one from definite conclusions. Below 
we will consider some group II NaNb0 3 - based solid solutions dielectric properties 
dependence on the concentration of the second and third components as well as on 
temperature and frequency. 

II EXPERIMENTAL RESULTS 

NaNb03-Gdi/3Nb03 solid solution crystals were grown by the flux method. The 
details of crystal preparation and characterization have been described elsewhere 
[7]. The ceramic samples of NaNb0 3 -based solid solutions were prepared by routine 
solid-state reaction route. The density of the obtained ceramics was about 90- 
95 % of theoretical one. For dielectric measurements Aquadag electrodes were 
deposited on the opposite faces of the as-grown crystals, while Ag paste was fired 
to the grinded disk surfaces of ceramic samples . The dielectric studies were carried 
out in the 1 kHz -1 MHz range in the course of both heating and cooling at a rate 
of 2 -3 C/min with the aid of the i?5083 and E7 — 12 capacitance bridges. 

Fig. 1 shows the typical evolution of the e'(T) dependencies with the compo- 
sition for the solid solution of group II. The similarity of the data obtained for 
ceramic samples and crystals shows that the peculiar properties of the group II 
solid solutions are not due to the immiscibility effect typical of many solid solu- 
tion ceramics. A strong frequency dispersion of both real (V) and imaginary (e") 
parts of complex dielectric permittivity observed at temperatures exceeding T m or 
(Fig. 2) is likely to be caused by a decrease of the conductance influence with 
increasing frequency. Though the e'(T) maxima of the group II NaNbC>3-based 
solid solution compositions with x > xq are smeared and a frequency dispersion of 
e' is observed, the shift of T m with frequency is usually much smaller than in the 
case of PbMg 1( / 3 Nb2/303-type relaxors (Fig. 2). Similar to other FE and AFE with 
DPT [1-3], the permittivity of NaNb0 3 - based solid solution compositions with 
x > Xq does not follow the Curie- Weiss law in a broad temperature range above 
T m (Fig. 2) while at higher temperatures the Curie- Weiss behavior is observed. It 



is well documented that in FE with DPT the extrapolated Curie- Weiss temperature 
of the FE phase transition, T C w, is necessarily higher than T rn due to a large con- 
tribution to the high-temperature permittivity from the regions that have a higher 
transition temperature [1-3]. In contrast to this, in all the group II NaNbOs-based 
solid solutions studied, the Tew values are much lower than T m (Fig.2). Figure 
3 shows the concentration dependencies of T m measured on heating for some solid 
solutions belonging to group II. It is interesting to note that for many solid solution 
systems of group II the linear extrapolations of the T m — x diagram portions from 
the x > x region intersect at the same point (x = 0, T ~ 150°C) corresponding 
to the minor e'(T) anomaly often observed in NaNb0 3 [8]. 
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FIGURE 1. Evolution of the e'(T) dependencies measured on heating and subsequent cool- 
ing with composition for (1-x) NaNbOa-xGd^NbOa single crystals (solid lines) and ceramics 
(dashed lines). Figures correspond to x values. 

A specific feature of the group II solid solutions is an anomalously large value 
of the thermal hysteresis AT/, of the e'(T) dependence typical of the composi- 
tions with x < xq. The ATh values increase with x and in some systems exceed 
100 K for compositions adjacent to x (Figs. 1, 4). The concentration depen- 
dence of AT h is substantially nonlinear (Fig. 4). When a small amount of FE 
third component is added to the given solid solution, ATh values usually increase 
[10-12]. At higher FE component content the ATh values decrease dramatically 
[11]. If the temperature is stabilized in the course of cooling within the thermal 
hysteresis range the temporal changes in e' do not exceed a few % during several 
hours [12]. Application of DC bias on compositions exhibiting giant ATh values 
does not change the e'(T) curve measured in the heating mode while in the course 
of subsequent cooling some lowering of the e'(T) maximim was observed [12,13]. 
Anomalously large values of AT h were observed both in the ceramics and single 
crystals and can serve as an experimental evidence of the fact that the given solid 
solution belongs to group II [7,9,11]. It is worth noting that in perfect NaNbOs sin- 
gle crystals the AT h value does not exceed 10 K , but increases substantially with 
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FIGURE 2. Temperature dependencies of e'(l'-3'), e" (l"-3") and 10 6 /e' (3) measured at 1 
kHz (l',l"), 10 kHz (2',2") and 100 kHz (3,3',3") for 0.88NaNbO 3 -0.12Gd 1/3 NbO 3 single crystal 




FIGURE 3. Concentration dependencies of the e'(T) maximum temperature T m measured on 
heating for some NaNb0 3 -based solid solutions belonging to the group II (1-5) as well as for 0.88 
NaNbi_ K Ta a; O 3 -0.12LiNbO 3 [16] system (6). The second components of the solid solutions are: 
l.BiFcOg; 2. SrCu 1/3 Nb 2/3 3 ; 3.CaTi0 3 [17]; 4. SrTi0 3 [23]; 5. NaTa0 3 [24]; 



the increase of the oxygen vacancy concentration [14]. In nominally stoichiometric 
NaNb0 3 ceramics the AT h values of 30 - 40 K are usually observed (Fig. 1). 
These larger values of AT/j in ceramic samples are likely to be attributed to a 
high point-defect concentration, e.g. due to evaporation of Na20 during sintering. 
Indeed, in Na20-deficient Nai_ x Nb0 3 _ :E /2 ceramics [15] the ATh values are sub- 
stantially larger than in the nominally stoichiometric ones. Both the oxygen- and 
Na20 -deficiency leads to an increase of the T m values in NaNbC>3 [14,15] and an 
almost linear correlation of AT h and T m is observed. 
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FIGURE 4. Concentration dependence of the s'(T) maximum temperature T m measured on 
heating (1) and on cooling (2) as well as the e'{T) temperature hysteresis AT/, values (3) for 
(l-x)NaNb0 3 -xGd 1 / 3 Nb0 3 crystals (open symbols) and ceramics (filled symbols) 
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FIGURE 5. Changes in the e'( T ) dependencies of 0.8NaNb0 3 -0.2Nao.5Bio.5Ti0 3 (1) and 
NaNbo.45Tao.55Os (2) solid solution ceramics caused by addition of 10 mol% of LiNb0 3 (curves 
la and 2a) or KNb0 3 (lb). Curve 2a is drawn using the data of Ref. [16]. 

It should be mentioned that a noticeable step in the T m (x) dependence of the 



NaNb0 3 -based solid solutions belonging to group II is observed only for the T rn val- 
ues measured in the heating mode. If T m is determined from the e'(T) dependence 
measured in the course of cooling, the T m (x) curve is rather smooth (Fig. 4). The 
e'(T) dependencies measured in the cooling mode are usually more diffused as 
compared to those measured upon heating (Fig.l). Moreover, just at the step in 
the T m (x) dependence the e'(T) curve obtained on cooling down has a relaxor- 
like diffuse maximum that can imply that the relaxor-type state had been already 
appeared in a metastable thermodynamic state occupied at cooling. The stability 
of this state arises at temperatures lower than the step in the T m (x) dependence. 
In the latter case (i.e. for the compositions with x > x ) the curves obtained in 
the cooling and heating modes practically coincide. 

As one can see from Figs. 1 and 3, the T m values of the known binary NaNb0 3 
-based relaxor-type solid solutions, belonging to group II, are usually well below 
the room temperature, and the maximal e' values, e' m , are much lower than in the 
Pb-containing relaxors. However, in the ternary NaNb0 3 -NaTa0 3 -LiNb0 3 solid 
solution system a dramatic increase in both the T m and e' m values with the 
LiNb0 3 content is observed [16] for the compositions from the x > Xq range 
(Figs. 3 and 5). Similar behavior is typical of other NaNb0 3 -based solid solutions 
(Fig. 5). 



Ill DISCUSSION 

When discussing the experimental data we will use a phenomenological Landau- 
type approach. For the sake of simplicity we do not follow the complex symmetry 
of the NaNb0 3 -based solid solutions but rather, in order to make a qualitative 
description, we consider a simplified case when there are two nonpolar (AFE) (Qi, 
and Q2) and one polar (FE) (P) order parameters. An important constraint on 
the Free energy expansion is that an AFE phase transition may exhibit a peak or 
step in the e'(T) only because of coupling between FE and AFE order parameters 
[18,19]. The impurities we consider are expansive that implies that they produce 
local random stress (07 ) leading to a local expansion of the lattice parameter and, 
hence, to local strain (e). The small observed dispersion of e'(T) at x > x requires 
the introduction of an additional parameter (d), which is the local mean square 
impurity- induced dipole (or multipole) moment. We assume that local expansion 
around a defect may lead to noncentrosymmetric relaxations of neighboring atoms 
yielding a net local moment, or a multipolar set of local moments. Averaged over 
the whole crystal, these relaxations do not generate a net macroscopic moment but 
local moments cause dispersion of e'(T). Extended dipole-dipole interactions are 
assumed to be sufficiently weak that they can be ignored. 

We expand the Free energy as a sum of x-independent (Fi) and rc-dependent 
(F 2 ) contributions, F = F 1 + xF 2 where 



F l = F + \a !(T)P 2 + i^P 4 - PP+ 

+\cx 21 {T)Q\ + \a 22 (T)Q\ - \p 2l Q\ + \(3 22 Ql + | 7 2iQ?+ 

+I/5mP 2 g? + |/9i22P 2 gi (1) 
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^2 = | c L e2 - °7 e - |aieP 2 + \a 2l eQ\ + \a 22 eQ\ - Arfe + \nd 

Here ctj = C\(T — T^), a 2 i = c 2 i(T — T$), and T^ - 1 are the bare Curie 

temperatures for the FE and for AFE phase transitions respectively, E is external 
electric field, e is the symmetrical part of the strain tensor. From the equilibrium 
condition with respect to e and d one has 

e = i ( ai + Xd + i fllP2 ~ l fl21 ^? ~ l fl 22^2) (2) 
d = Ae/re 

It implies that d multiplied by A plays the role of the local stress and hence it 
effectively enlarges the local stress produced by impurities, d proves to be propor- 
tional to the strain e. From (1) and (2) one can easily find that e (as well as d, 
see [20]) is enhanced due to mutual coupling between e and d: 

e = ^ ^ (a t + \a x P 2 - \a 21 Q\ - \a 22 Q$) (3) 

By using (3) one can rewrite Free energy in the simple form containing expansions 
only with respect to P and Q: 

F = F[ ) + \A X (x, T) P 2 + \B X (x) P A -EP + \A 21 {x, T) Q\+ 

\A 22 (x, T) Q\ - \B 2l (x) Qf + |r 2 (x) Q\ + \B 22 (x) Qi+ (4) 

\B l2l (x, T) P 2 Ql + \B l22 (x, T) P 2 Q 2 2 + ... 

where Fq = P — xaf/2q and coefficients A and B depend on concentration, e.g. 

B 21 (x) = p 21 + xal 1 /2<; 

B 121 (x) = fczi + xaia 2 i/2q (5) 
^122 (x) = P122 + xa x a 22 /2q 

Here q = (cl — A 2 //t) 2 /cl- We have obtained that the coupling of the order pa- 
rameter Qi with the strain e produced by the impurities leads to the dependence 
of the coefficients in the Landau expansion on the concentration x (see also [21]). 

Composition dependence in coefficients (5) implies that the Curie temperature 
also depends on x: 

Ti = T 1 (0) + xaiaxjcxq , ■ 

rp rr(0) I W 

T 2i = T 2i ~ xaia 2i /c 2i q 



From these expressions it is seen that if a\ > and a 2 > then Ti increases 
with x and T 2 decreases with x linearly. The latter result corresponds to the 



experimental data: all group II NaNb0 3 -based solid solutions studied exhibit a 
decrease in T rn with x. The assumed signs of a\ and a 2 are natural as the 
electrostrictive constant a± is always positive (if in the F one puts the sign "-" in 
front of it) and a 2 should also be positive (if one puts the sign "+" in front of it) 
because compression of the lattice (e < 0) typically promotes AFE order while the 
expansion (e > 0) promotes FE order. 

From expansion (4) one can easily obtain: the critical temperature of the AFE 
phase transition, Tau the temperature hysteresis width, AT h ; and the jump of the 
AFE order parameter, AQi. 

Tax = T 21 + 3B 2 J16c 2ll21 - B 121 P 2 /c 21 
Ta2 = T 22 — B\ 22 P 2 1 c 22 
AT h = B 2 21 /4c 2ll21 

AQi = ^21/2721 

From (6) it is clear that T 2 « is a linear function of x, which implies that T A i{x) is 
also linear at small x. This corresponds to the experimentally observed linear 
decrease of T m {x) with x both above and below x as well as explains why 
linear extrapolations to the point x = for various second components, intersect 
at one point (Fig. 3). These points correspond to the bare Curie temperatures 
T 2 ®\ The slopes of the linear dependencies above and below x are different 
as are the electrostriction constants (a 21 > a 22 ). Thus: Curie lines T A i{x) and 
Ta2(x) intersect at x — xq] Q\ is nonzero for T < Tax (x < %o)', Qi is nonzero 
for T < Ta2 (xo < x ). The point x = Xq can be found from the equality: 
T A i(x) = T A2 (x). 

At the intersection of curves T A i(x) and T A2 (x), the first order phase transition 
to the AFE phase, with order parameter Q X) intersects the second order phase tran- 
sition to the AFE phase, with order parameter Q 2 , at a critical endpoint [22]. The 
principle experimental evidence supporting this interpretation is that the thermal 
hysteresis vanishes abruptly at x (Fig.4). Note that we assume B 2 \ > and 
B 22 > but the sign "-" stands in (4 ) in front of B 2 \ and "+" precedes B 22 . 
The alternative assumption, that there is a tricritical point at (xo, Tai{x)\ would 
imply a gradual disappearance of the thermal hysteresis as x approaches x from 
below. 

Strong rc-dependence of the thermal hysteresis at x < x is unusual. From (5) 
it follows that B 21 increases with x regardless of the sign of the electrostriction 
constant. At small x, this leads to a linear increase in AT/,(x) while at larger x, 
ATh(x) increases quadratically, in excellent agreement with experiment (Fig.4). 

From (7) it follows that the jump, AQi, in AFE order parameter Qi should 
depend on the concentration x. Because B 2 increases linearly with x, so must 
AQx. Experimentally, it is observed that as x increases, so does the difference 
between the magnitudes of e' m measured on cooling and on heating (Fig.l), which 
implies good agreement between theory and experiment. 
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Our data are insufficient to identify the crystal structure of the AFE phase asso- 
ciated with Q 2 = finite (xo < x) but some properties are clear: e'(T) is relaxor- like 
with a diffuse maximum and dispersion below T m . Some properties however, are 
not typically relaxor- like: the dielectric permittivity magnitudes are much lower 
than those in typical relaxors such as PMN; the e\T) fits the Curie- Weiss law 
well above T m , but the extrapolated Curie- Weiss temperature, T C w, is significantly 
below T m rather than above as in typical relaxors; in spite of the dramatic diffusion 
of e'(T), the frequency dependence of e' and T m is much weaker than in typical 
relaxors. The comparatively small magnitude of e' m can be connected with the 
absence of the large contribution of lead to the dielectric permittivity. 

IV SUMMARY 

From the discussion above it follows that one of the main reasons for the appear- 
ance of the relaxor-like properties in the NaNb03-based solid solutions belonging 
to group II can be the appearance of the (random) local strain (with nonvanishing 
average magnitude) stemmed from the impurities. Due to the electrostriction effect 
this results in the decrease of the AFE critical temperature, increase of the thermal 
hysteresis width with x, and, finally, with the appearance of a new diffuse phase 
transition with weak dispersion of the dielectric permittivity. Large diffusion of the 
e'(T) maxima in some solid solution compositions with x > x and the possibility 
of shifting T m to the room temperature range in conjunction with relatively weak 
frequency dependence of both T m and e' m may be of certain interest for appli- 
cations. The same is true for solid solution compositions with x < xq exhibiting 
giant e'(T) thermal hysteresis values. 
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